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Abstract 

A procedure to characterize chaotic dynamical systems with concepts 
of complex networks is pursued, in which a dynamical system is mapped 
onto a network. The nodes represent the regions of space visited by the 
system, while edges represent the transitions between these regions. Pa- 
rameters used to quantify the properties of complex networks, including 
those related to higher order neighborhoods, are used in the analysis. The 
methodology is tested for the logistic map, focusing the onset of chaos and 
chaotic regimes. It is found that the corresponding networks show distinct 
features, which are associated to the particular type of dynamics that have 
generated them. 

89.75.Fb: Structures and organization in complex systems 
89. 75. He: Networks and genealogical trees 
02. 10. Ox: Combinatorics; graph theory 



1 Introduction 

Chaotic dynamical systems are characterized by several measures that quantify 
how irregular (despite deterministic) the trajectories are. The set of Lyapunov 
exponents [T] provides a measure of the dependence on the initial conditions 
of their trajectories, while information theory may be used to characterize a 
dynamical system in terms of the production of entropy. Actually, a dynam- 
ical system with a chaotic behavior is regarded as a realization of Shannon's 
concept of an ergodic information source [5], as the Kolmogorov-Sinai entropy 
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is equated [3] to the sum of the positive Lyapunov exponents. Further mea- 
sures to describe a chaotic system include fractal dimensions and singularity 
spectra [H0. This formalism applies only where the system has, at least, one 
positive Lyapunov exponent. However, several authors have studied situations 
in which the largest Lyapunov exponent vanishes, but the distance between 
neighbouring trajectories increases in time according to a power law [5J [7J [5] [H] ■ 
These situations are usually found at the onset of chaos, where an infinitesimal 
change of a control parameter drives the system into either a regular or a chaotic 
regime. These investigations uncovered some of its features related to the sen- 
sitivity to the initial conditions [5J, entropy production per unit time [TO], 
multifractal geometry of the attractor [TTJ, relaxation to the system attractor 
[T2] and multifractal dynamics at the onset of chaos [13] . 

Recently, the investigation of complex networks has set up a new framework 
for the analysis of systems with a large number of degrees of freedom. Within 
it, one has access to the properties of the topological structure underlying the 
mutual interactions among the system constituents. This approach has been 
applied to a large variety of actual systems, ranging from social interactions, 
biological nature, information and electrical distribution [14l [15] . 

In this work, we define a mapping of dynamical systems onto a network. 
This way, the network properties can be used to display new features for the 
characterization of the system's trajectory. The network nodes correspond to 
coarse-grained regions (cells) of the phase space visited by the trajectory. Two 
nodes r and s are linked if, during the time evolution, the trajectory jumps from 
cell r to cell s. Although this naturally offers a construction procedure for a 
directed network, we consider here undirected networks. In this approach, we 
are able to find novel geometric and topological properties of the phase space 
through the measures that have been recently developed to characterize com- 
plex networks. The dynamical system, defined on the time domain, is mapped 
into a node domain, which represents regions of the phase space visited by the 
trajectory. As it is based on the division of phase space into boxes, this pro- 
cess follows similar construction procedures considered in the evaluation of the 
fractal dimension of attractors and the Kolmogorov-Sinai entropy. 

It worths mentioning that some previous works have put together ideas of 
dynamical systems and complex networks. As we shall see, our approach is 
rather distinct, as the previous contributions consider synchronization under the 
assumption of certain regularity in the connection topology |16[ I17j . The use 
of a network to represent the phase space evolution of discrete time dynamical 
systems has been suggested in [T5] , 

We restrict ourselves to the analysis of the quadratic map [19] described by 

x t+1 = \-ax\ (t = 0,1,2,3- ••) (1) 

where x t £ [—1,1] and a £ [0, 2]. Eq. (JTJ) leads to a variety of distinct dynamical 
situations, the properties of which are expected to be manifested in the networks 
they originate. In particular, we explore those regimes at the onset of chaos that 
proceed through a bifurcation cascade as well as those close to an intermittency 
transition and also the full developed chaos. 
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2 Network characterization 



An indirected network R is denned only by the number N of nodes and links 
L, represented by an assembly of unordered pairs Sn(r, s), r, s < N, indicating 
which pairs of nodes are directly connected. This information provides a full 
description of the network, leading to the computation of the average number 
of links per node (fc), the average clustering coefficient C, the mean minimal 
distance among the nodes (d), the diameter D, the probability distribution of 
nodes with k links p(k). Other measures, like the assortativity degree [20] and 
the distribution of individual node clustering coefficients C{k) with respect to 
its degree k have also been introduced, but they are not explored here. 

R can be described by its adjacency matrix M(R). This is not the most con- 
cise representation of a network, but it opens the possibility to the evaluation of 
its spectral properties and, as recently indicated, the higher order neighborhoods 
Ri, £ — 1,2, D [5T]. This is done in a straightforward way, by means of the 
set of matrices {Me}, so defined that (Me) rjS = 1 only if the shortest distance 
along the network between the nodes r and s is i. Otherwise, (Mi) r:S — 0. Al- 
though the whole information on the network is entailed in Sn(r, s) or in M(R), 
each Mi condenses information on R that is extracted from M(R) within the 
quoted framework. This formalism is also consistent with the recently proposed 
procedure to evaluate the fractal dimension of the network (Lf,r }22j . as it nat- 
urally leads to the set {Ng} required for this evaluation. Here, each Ne counts 
the number of pairs of nodes which are £ steps apart. 

Within this framework, we consider that each node is the only zeroth order 
neighbor of itself, and define Mq = I, where / indicates the identity matrix. 
Also, we assume that M\ — M. Since the matrix elements of M(R) take only 
the values and 1, the other matrices Mi of the set are recursively evaluated 
with the help of Boolean operations [3T] . We shall make further use of a matrix 
that condenses all information in {Me(R)}. As just discussed, given any two 
nodes r and s, it is clear that {Mg) rs — 1 for just one value of £. So, if we define 
a matrix 

M= y £ f jM J , (2) 

3=0 

it directly informs how many steps apart any two nodes in the network are. Also, 
it is possible to use the information in M to graphically display the structure 
of a network with the help of color or grey scale plots. 

It should also be mentioned that this framework opens the door for a finer 
characterization of the network, if we consider each Ri as an independent net- 
work. Thus, several of the properties quoted above used to characterize R can 
be also evaluated for the evaluation of Ri . This is explored in the next section 
specifically for the the degree distribution and clustering coefficient (k)e and 
C{£). 



3 



3 An algorithm to map a dynamical system into 
a network 



To map a dynamical system into a network R, we use the framework of the 
algorithm introduced in 23J, to efficiently evaluate the generalized fractal di- 
mensions of fractal structures by the box counting method. 

Let us consider a dynamical system with m variables. Without loss of gen- 
erality, one can consider a set of points Z C 5ft m consisting of the vectors z(i), 
!=l,-",T,r>l, which represent the coordinates of the dynamical system. 
The components of these vectors, Zs (i), 5 = 1, • ■ • , m, are assumed to belong to 
the interval [0, 1 ). We define a graining in phase space by dividing each phase 
space axis into W equally sized disjoint intervals, so that the whole phase is 
spanned by a set of W m boxes. This represents also the maximal possible num- 
ber of nodes in a network, as e.g., in the case of an ergodic system. Of course, 
the choice of W defines the graining, and the size of the region represented 
by a node. In the next section we investigate the effect of W on the obtained 
networks. 

Based on |23j . each point z(i) of the trajectory is mapped into a node of R 
according to 



where floor(x) is a function that evaluates the largest integer less than x. Actu- 
ally, this is just a simple way to divide the region [0, l) m in equal parts. Thus, 
the nodes of the so constructed network represent a box in the coarse grained 
phase space of the system. After the mapping is complete, the boxes that were 
not visited by the trajectory are eliminated from the network, as they constitute 
nodes with zero degree (k = 0), which do not provide any useful information on 
the dynamical system. The edges are built as described in the following proce- 
dure. Let z(i) and z(i + 1) be two consecutive points of the dynamical system. 
Suppose that these points were previously mapped into the nodes n(r) and n(s), 
where < r, s < W m . Thus, one assumes that there is an edge connecting n{r) 
and n(s). Here, it is considered that there is only one edge linking n(r) to n(s) 
if r =/= s, i.e. self-links are not computed. 

This procedure can lead to directed and weighted networks; however, in this 
work, we focus on the most simple situation of undirected and un-weightcd 
networks, once our main purpose here is to address the problem, and show that 
we can extract useful information from it. 



We concentrate our investigation on three distinct regions of a, a = a c — 
1.40115518909..., a G [1.749, 1.75) and a = 2, which correspond, respectively, to 
the first period doubling transition, the region close to the tangent bifurcation 
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4 Results 



4 



to the period-three window, and the full developed chaotic state. Representa- 
tive networks to the distinct chaotic attractors are generated for distinct values 
of the graining W. We only consider trajectories that start on the attractor, in 
order to avoid spurious nodes (visited only once) that depend on initial condi- 
tions. For a fixed value W, the network grows as the trajectory evolves in the 
the phase space with respect to the number of iteration steps t. There is no a 
priori criterion to decide the time tp after which the network is complete. In this 
work we have followed how N and L increase with t, for a given W. We define 

as the smallest value of t for which N(tp) — N(2tp) and L[tp) — L(2tp). 

First let us discuss the effect of W on N and L. For the purpose of presenting 
a full neighborhood analysis of the networks, we have selected here values of W 
that lead to the maximal number m 10000 nodes in the network. Such choice for 
W clearly depends on a. Indeed, due to the strategy adopted for the construction 
of the networks, N grows with W according to a power law mediated by the 
fractal dimension of the attractor dp,A- This is shown in Fig. [TK, where we draw 
points for a = a c , a c +10~ 3 , 1.749999 and 2. For a c we measure the slope 0.54 . . ., 
which agrees with the known value of (1f,a of the period doubling attractor. For 
all other cases, the slope is 1 within 2% accuracy, even for a = a c + 10 -3 , which 
lies already in the chaotic regime. This is in accordance with the fact that dp t A 
changes in a discontinuous way at a — a c . When a corresponds to a periodic 
solution, the network becomes finite, so that N and L do not depend on W, 
provided this parameter is large enough. 

Although we are primarily interested in the properties of the complete net- 
work, we can also follow the dependency of N and L with respect to t, for a fixed 
value of W. Assuming N ~ L z , this defines a dynamical exponent z in the early 
stages of the network evolution. The analyzed data indicate that z ~ 1 for all 
values of a. Nevertheless, we find that, in the immediate chaotic neighborhood 
of a = 1.75, the laminar phases in the intermittent regime traps the trajectory 
for long intervals, demanding large time of integration to complete the network. 

Now we discuss the network properties, following the methodology and pa- 
rameters indicated before. When appropriate, we extend our discussion to prop- 
erties of higher order neighborhoods in the network. We find distinct network 
structures when we consider the chaotic regime or the onset of chaos. Regarding 
the mean minimal distance (d) and diameter D, we find that, for the chaotic 
regime, they grow with respect to W (and N), in a logarithmic way, similarly 
to small- world networks 124:. If we write (d) = alog 10 W, we find a ~ 2.4 and 
3.7, respectively, for a = 2 and a S [1.749, 1.749999], as illustrated in the inset 
of Fig. [TJd. As D only assumes integer values, it is possible to expect a similar 
behavior only in an approximate way, e.g., equally sized steps in D x log 10 W 
plots. Thus, assuming D — /31og 10 W, we find f3 ~ 4 for a = 2 with very good 
accuracy. In the interval [1.749, 1.749999], we noticed the presence of fluctua- 
tions in the size of the steps, which increases when we approach the threshold 
a = 1.75. The results for a = a c have a completely distinct behavior: (d) and 
D increase as power laws with respect to W, as illustrated in the main panel 
of Fig. Ql>. For a — a c + 10~ 3 , the same dependence prevails. The exponents 
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Figure 1: (a) Dependence of N with respect to W for a = a c (circles), a = 
a c + 10~ 3 (diamonds), a = 1.749999 (down triangles), and a — 2 (up triangles). 
The the same convention is used in all other figures. The distinct slopes indicate 
the values of cIf,a- (b) Power law dependence of (d) (hollow symbols) and D 
(solid symbols) with respect to N for a = a c and a = a c + 10 -3 . The inset 
shows logarithmic dependence among the same quantities when a = 1.749999 
and a = 2. 
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Figure 2: Degree distribution of nodes at a c (circles) and in the a — 2 chaotic 
regime, (up triangles). Behavior close to tangent bifurcation, a = 1.75 — 10~ 4 , 
(not shown) is quite similar to a = 2. 

obtained for (d) and D are, respectively, 0.67 and 0.73 for a c , and 0.35 and 
0.38 for a = a c + 10~ 3 . (d) and D increase smoothly for a = a c + 10 -3 but 
at a = a c , the growth of the two quantities present discontinuities and steps. 
Results in Fig.[T)3 indicate that, unlike the properties of the attractor, reflected 
by the rash change in the value of dp. a, the network properties change slowly 
when the chaotic regime is reached. 

We have evaluated the distribution p(k) vs. k in all different regimes. We 
can see in Fig. [2] that, for a = a c , k does not reach large values (k max ~ 30), 
so that it is not possible to identify a power law decay in this range. For a = 2 
(and also a = 1.749999) nodes with larger values of k can be found, but p(k) 
does not follow either a power law. 

Such distinct behavior is also explicit when we analyze the fractal dimension 
of the network dp,R [22- Fig. [3] shows that the a = a c networks have a well 
defined scaling behavior, which extends over more than two decades, in a quite 
precise way. On the other hand, a fractal dimension for the networks in the 
chaotic regime is not evident. First, the small values of D reduces the region of 
possible scaling behavior. Then, we clearly observe deviations to the expected 
power-law regime. 

Regarding the clustering coefficient, we obtain C = for a = a c , indicating 
the complete absence of triangles in the network. For a = 1.749999 and a = 2, 
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Figure 3: Clear power law behavior for N(£) x I when a — a c , with dpn = 1.47. 
Finite size effects blur this dependence when £ ~ D. In the chaotic regime, for 
a = 2, <If,r can hardly be evaluated. For a = a c + 10~ 3 , the points illustrate a 
slow transition between the two regimes. 
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we find that C decays with iV according to a power law with exponent w 0.95, 
what is not so close to the value 0.75 observed for the Albert-Barabasi scale- free 
network [T2]. However, its small values indicate that the network has only a 
small number of triangles. 

Other features of the network may be drawn if we consider the clustering 
coefficient of higher order neighborhoods [5T] . To obtain a clearer picture of this 
analysis consider, for instance, the regular Cayley tree with 1 = 2 network for 
the complete regular three neighbors Cayley tree. The 1 = 1 network is formed 
by triangles, much like the Husimi cactuses and, correspondingly, a large value 
for C{£ = 2). The following odd and even numbered neighborhoods are charac- 
terized, respectively, by values ofC(£) = and C(£) > 0, whereby the values of 
C (I) for a subset of even neighborhoods decrease in a monotonic way. A similar 
situation is found for the networks investigated here. In Fig. [4^, we summarize 
the sequence of C(£) for the three situations under investigation. We see that, 
for a — 2 and a — 1.4999, the oscillatory behavior lasts only until £ = 5 and 
10, respectively. Also we note that the odd numbered C(£) increases until it 
reaches values as high as those of the even numbered neighborhoods. On the 
other hand, the a = a c network has C(t) = for all odd numbered neighbor- 
hood. The inset shows that the C(£=even number) decays with I according to 
a power law, with exponent a ~ 1. 

We have also analyzed the average degree (k)e as a function of the neigh- 
borhood £. Here again we find that the behavior of the chaotic regime and the 
onset of chaos have distinct features, as exhibited in the Fig. [4)d. 

Finally, we use the information in M to obtain images, in color (or gray) 
scale, of the network neighborhood structure. They provide a vivid and easy 
visualization of their distinct properties of the attractor. For a = 2, the first 
order neighborhood is distributed along the parabola described by the r.h.s. of 
(JTJ (see Fig.[5Ji). It shows how the second and higher order neighborhood evolve 
according to the higher order iterates of the quadratic map. However, mixing 
and finite size effects stemming from a finite graining smears the higher order 
iterates. The situation is distinct for a c , when the attractor is a dp a — 0.54... 
dust spread out in the [—1,1] interval. Only boxes containing part of the dust 
remain in the network, so that contiguous numbered nodes are not actually 
neighbors in phase space. The resulting image (Fig. [5Jd) displays a fine and 
intertwined tessitura that reflects a very peculiar behavior of the trajectories at 
the onset of chaos [3 EE]. 

5 Conclusion 

In this work we explore the idea of mapping chaotic systems onto complex 
networks. Networks are constructed according to a well defined methodology, 
and results for the logistic map indicate how their properties can be associated to 
those of the attractor in phase space. Trajectories in distinct dynamical regimes 
are explored in order to show how the major differences in phase space are 
reflected in the networks. The networks show several features of small- world and 
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Figure 4: (a) Behavior of C(£) with respect to £ for chaotic regime and at the 
onset of chaos. The inset shows a peculiar power law behavior C(£) at the onset 
of chaos, (b) Behavior of (k)i x £. The onset of chaos exhibits a slow increasing 
value of (k)i over a large interval of £, interrupted by finite size effects already 
present in Fig. [3] In opposition to this picture, chaotic regime shows a sharp 
increase of (k)e- 
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Figure 5: Color scale plots (gray scale in paper version) of M for a = 2 (a) and 
a = a c (b). Scale ranges from black and blue (£ = and £ = 1) to red (£ = D). 
Number of levels in (a) is much smallpj than in (b). Neighborhood structure 
changes abruptly in the distinct regimes. 



scale-free networks, but they do not fully match with those networks generated 
by the specific algorithms described in [24l [25] . The analysis of the networks in 
the neighborhood of a c reveals that the N xW dependence, measured by dp t A, 
has a sharp transition at the onset of chaos. So, the distinct character of the 
trajectories in phase space is indeed reflected in the network. However, with 
exception of C, the results for the other indices {(d), D and p(k)) change in a 
much smoother way with respect to changes in the parameter a. 
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